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Summary 

This report presents a method for estimating the 
crack-extension resistance curve (R-curve) from 
residual-strength (maximum load against initial crack 
length) data for precracked fracture specimens. The 
method allows additional information to be inferred 
from simple test results, and that information can be 
used to estimate the failure loads of more 
complicated structures of the same material and 
thickness. 

The analytical basis for the method was developed 
by the author in a previous report. Numerical 
differentiation of the residual-strength data is re- 
quired, and the problems that it may present are 
discussed. The method is demonstrated and verified 
by using various types of residual-strength data from 
the literature. 


Introduction 

This report presents a method for estimating the 
crack-extension resistance curve (R-curve) from 
residual-strength (maximum load against original 
crack length) data for precracked fracture specimens. 
Such data may be found in the literature or may exist 
in company files. The method allows additional 
information to be inferred from simple test results, 
and that information can be used to estimate the 
failure loads of more complicated structures (such as 
reinforced panels) of the same material and 
thickness. 

The progressive development of the R-curve 
concept has been reviewed in reference 1 . The 
concept postulates that, for a given material and 
thickness, there is a unique relationship between the 
amount of stable crack growth under rising load and 
the crack-tip stress intensity factor. The relationship 
is called the crack-extension resistance curve, or 
R-curve, and represents the response of the material 
near the crack tip to externally imposed loading. If 
the R-curve is known, both failure load and critical 
crack length can be predicted (as functions of initial 
crack length) for any specimen or structural 
configuration for which an appropriate stress 
intensity analysis is available. For example, refer- 
ence 2 illustrates the use of the R-curve to predict the 
failure load for a strap-reinforced panel. Current 
methods for determining the R-curve (ref. 3) require 
that the specimen be instrumented to measure crack 


opening displacement, and the R-curve is derived 
from the load-displacement record. 

The analytical foundation for the present report is 
developed in reference 4. That report examines 
several semiempirical fracture analyses that attempt 
to relate maximum load to initial crack length in 
terms of one or two empirical parameters. It showed 
that an R-curve can be developed that is 
mathematically equivalent to any one of the semi- 
empirical analyses. In the present report the R-curve 
is developed directly from residual-strength data 
without recourse to any semiempirical analysis. 

This report first reviews some characteristics of the 
R-curve concept when applied to practical test 
specimens. Next the method for generating an esti- 
mated R-curve is derived. Numerical differentiation 
of the residual-strength data is required, and the 
problems that it may present are discussed. Finally 
R-curves are estimated for several materials by using 
data from the literature. 


Symbols 

a length of single-tip crack or half-length of 

double-tip crack 

E' effective modulus, E for plane stress or 

E/(l — v 2 ) for plane strain, where E is 
Young’s modulus and v is Poisson’s ratio 

G a strain energy release rate 

G c fracture toughness, G A or Gr at instability 
condition 

Gr crack extension resistance 

Kj opening-mode stress intensity factor 

n number of crack tips (one or two) 

W specimen width 

Y stress intensity calibration factor, Kj/a^fa 

a sensitivity factor (eq. (1)) 

7 sensitivity factor (eq. (5)) 

A effective crack extension (sum of physical 

crack extension plus a plastic zone 
correction) 

X relative crack length, na/W 

a stress normal to crack 

o u ultimate tensile strength 

o ys yield strength 


Subscripts: 

c at critical or instability condition 

0 initial value (before loading) 

R-Curve Concept 

The R-curve concept is illustrated schematically in 
figure 1 for an infinite body containing a crack whose 
original length is 2oq . The strain energy release rate is 
given by 

o 2 ira 

GA =~E r 

and represents the driving force (per unit thickness) 
tending to cause crack propagation. The material’s 
resistance to crack propagation Gr is a function of 
crack extension A. In R-curve analysis the subscripts 
A, R , and c are customarily used to denote applied 
and resisting forces and critical values, respectively. 
The R-curve is located with its origin at a=ao- As 
stress normal to the crack is applied and increased to 
90 percent of the subsequent critical stress in figure 1, 
the crack must extend only a small distance to 
develop a large resistance. At this point the crack- 
extension resistance equals the driving force, and the 
crack is stable. As the stress is increased. 



Crack length, a 

Figure L - Schematic representation of R-curve in- 
stability concept for an infinite body. 


progressively larger amounts of crack extension are 
required to resist the crack driving force. Finally at 
the critical stress a c the driving-force curve and the 
R-curve are tangent. Beyond the point of tangency 
the driving force increases faster with crack length 
than does the material’s resistance. This instability 
condition represents the failure of the body. The 
point of tangency defines the fracture toughness G c 
and the critical crack length 2 a c . Since the driving- 
force curve for an infinite body is a straight line, it 
should be apparent that both the fracture toughness 
G c and the amount of crack extension at instability 
A c increase with increasing original crack length 2oq . 
If the R-curve exhibits a plateau, G c and A c may 
asymptotically approach limit values. 

Instability Calculations 

In simple finite bodies and test specimens the 
presence of stress-free boundaries results in an 
additional increase in the crack driving force as the 
crack extends toward a boundary. Thus the slope of 
the driving-force curve increases continuously with 
increasing crack length. The instability condition for 
a typical finite-width specimen is shown in figure 2 
and is determined as follows. For a given specimen 
configuration and loading the dimensionless stress 
intensity calibration factor is defined as 

Y~ =f cn (\) 

cNa 

where Kj is the opening -mode elastic stress intensity 
factor, X^nafWh the relative crack length, and n is 
the number of crack tips (one or two). If we define a 
dimensionless sensitivity factor as 



Figure 2. - R-curve instability concept for a finite body. 
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X dY 
a ~ Y dk 


equation (2) by using the appropriate value of Y c . In 
W this report the expression 


then (for constant stress) the crack driving-force y=(?r secant — )* /2 
curve and its slope are 2 


E'G a = Y 1 o 2 a 


is used for the center-crack specimen (ref. 5). 


E'^d. =Y 2 a 2 d+2a) 

da 

For convenience the crack-extension resistance curve 
and its slope are written here as 


g(A)=E'GR 


g'(A)=E 


, dG R 
dA 


At the instability point, G A = Gr and 
dG A /da=dG R /dA (fig. 2). If g(A) and g'(A) are 
mathematically describable, the instability point is 
determined by the simultaneous solution of the 
following equations: 


E'G c = Yy c {a 0 +A c )=g{A c ) 


( 2 ) 


Instability in Finite Bodies 

For a crack in an infinite body both G c and A c 
increase with increasing initial crack length. In a 
finite body as the initial crack length is increased 
from zero, G c and A c increase at first. However, 
because dY/cf\ continually increases with X, both G c 
and A c reach maximum values that depend on the 
specimen width W and the forms of both the driving- 
force curve and the R-curve. As the initial crack 
length is increased still further, both G c and A c begin 
to decrease. This behavior is shown schematically in 
figure 3, where instability curves are shown for a 
wide range of initial crack lengths. The locus of all 
instability points is shown by the dashed line. 

It is useful to define a sensitivity factor 


1 + lot __ n 1 dG A 
X ~WG^ 'da 


(5) 



= y^(l+2a c ) = g'(A c ) 

<*c 


which is shown for the center-crack specimen in 
(3) figure 4. For comparison, y is also shown for the 
standard compact (tension) specimen. This was 
derived by using the expression 


The coefficients Y and a are usually expressed as 
trigonometric or polynomial functions of the relative 
crack length X. As a result a closed-form 
simultaneous solution is seldom possible, and 
numerical methods must be used to solve for A c . This 
can be done as follows: Dividing equation (3) by 
equation (2) and rearranging terms give 


_ g'(A c ) l+2* c 

g(A c ) a 0 + A 


If the functions g(A) and g'(A) and the appropriate 
equation for a are substituted into equation (4), then, 
for prescribed values of ao and W, A c is the least 
positive root of equation (4). This root can be found 
by any of several numerical methods. Next, the 
function g(A) is evaluated at A c to calculate g(A c ). 
Finally the fracture stress o c is obtained from 


Y= X- 1/2 (l - X)- 3/2 (2 4- X)(0.886 + 4.64 X 

- 13.32 X 2 + 14.72X 3 - 5.6 X 4 ) 



Figure 3. - R-curve instability for a wide range of initial 
crack lengths. 
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Figure 4 - Sensitivity factor for center-crack and com- 
pact specimens. 


from reference 6. Note that the function X has a 
minimum value. Then for a given specimen width the 
first term in equation (4) must also have a minimum 
value. In graphical terms this means that if, at a given 
point, the R-curve’s slope is too low or its magnitude 
too high or both, a driving-force ( G a) curve cannot 
be both coincident and tangent at that point. Thus 
for a given specimen width the value of X at 7 m j n 
defines the critical crack length at which both crack 
extension A c and fracture toughness G c or K c are 
maximum. Note that both 7 m j n and the relative crack 
length at which it occurs are different for different 
specimen types. The fact that the function 7 has a 
minimum has additional significance that is evident 
later in this report. 

Comments 

Although the subject has received scant attention 
in the literature, it is possible to estimate the R-curve 
from residual-strength data by using purely graphical 
methods. The method is shown in figure 5. To 
simplify the illustration, assume that we know the 
residual strengths of four specimens having the same 
initial crack length but different widths. Construct a 
driving-force ( Ga ) curve for each specimen by using 
residual strength in the equation following equa- 
tion (1). Then draw the estimated R-curve from the 
point (a = ao, Ga = 0) in such a way that it is tangent 



Figure 5. - Example of graphical method for estimating R- 
curve from residual-strength data. (W 4 < W 3 < W 2 < Wj = 
00 ; 04<a3<02<o^) 


at some point to each of the driving-force curves. If 
initial crack lengths are not constant, it is convenient 
to let the x-axis be given by crack extension A. Then 
each driving-force curve begins at a distance ao to the 
left of the origin. Although the procedure appears 
simple, in practice it is tedious and subjective at best. 
If significant data scatter exists, one’s skill with the 
French curve will be severely tested. 

In concluding this section it should be emphasized 
that, in this report, A is the effective crack extension. 
It is the sum of the physical crack extension plus an 
adjustment to account for the effect of crack-tip 
plasticity. The nature of the plasticity adjustment, 
although technically significant, has no influence on 
the analyses and conclusions in this report. 


Analysis 

The analytical foundation for this report was 
developed in reference 4. However, the common 
elements of the derivation are repeated here for con- 
venience. In the preceding section it was shown that, 
if a mathematical formulation for the R-curve is 
available, fracture stress can be predicted as a 
function of original crack length. In this section it is 
shown that, if proper estimates of residual strength 
and its first derivative can be made, the R-curve can 
be estimated. 

The first step requires that we differentiate 
equation (2) with respect to A c > and that operation 
requires some prior consideration. For a given spec- 
imen type there are many combinations of specimen 
width and initial crack length that will result in 
instability at the same point on the R-curve, and each 
combination has an associated fracture stress. To 
cause a small change dA c in the instability point, 
there must be a change in the specimen width or the 
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initial crack length or both, and there will be a 
change in the fracture stress. Thus, when 
differentiating equation (2), the terms W, o , and ao 
are treated as variables. 

Differentiating both sides of equation (2) with 
respect to A c yields 


dA c 


E’G c = Y 2 c a 2 c 2a c 


v( da ° _i_1 _ g 0 +A c dW \ 

^dA c W dA c > 


+ Y*(a o + A c ) -jl- (pi) 




and substituting equation (3) results in 




-2a, 


a 0 + A c dW 
W dA, 


Then equation (6) becomes 
0 = ^ [ (1 + 2o c V(a 0 ) + (ao + A c ]f' (fl 0 )] 

and, since da^/dA c ^0, we have 


0= 1+2< *c | /^gQ ) 
«o+ A c A«o) 


(7a) 


Case II: ao/W = Constant [dW = (W/ao)dao]. — 
Now assume that there is a function h such that 


h(a 0 )=o 2 c 


'(n«\= tf ( g c) 


A'(«o)« 


ofao 


a 0 

for — = Constant 

w 


J 


Now equation (6) becomes 


°=^ ( J ~2a c ^)^(a 0 ) + (ao +A c )/!'(a 0 )] 


and since da§/dA c ^0, we have 


+ ( a o + A c) (o’) 


(6) 


At this point it is helpful to reduce the problem to 
one of a single independent variable by prescribing 
the manner in which ao and W may vary. Three cases 
are considered. 

Case I: W = Constant (dW = 0). — Assume that 
there is a function / such that we can define 


f(ao) = al 


-\ 


for W~ Constant 


/'<«<>)- 


d ( ff c) 

dao 


Q = 1 — 2a c (A c /a 0 ) ^(gp) 


(7b) 


Case ///.* ao —Constant (dao =0). — This case is of 
limited usefulness but is included for completeness. 
Assume that there is a function j such that 


j(W)^a 2 c 1 


J'W 


_ d {°V) 


dW 


for ao = Constant 



Now equation (6) becomes 


0=^(a 0 + A c) 


2a c 

~W 


jm+j'w 


and since d\V/dA c 9±0 and (ao + A c )^0, 


2a c J'W 


(7c) 


In reference 4 the functions /, h, or y and their 
derivatives were obtained from analytical 
formulations presented in various semiempirical 
fracture analyses. However, if sufficient and suitable 
experimental data are available, these functions and 
derivatives can be obtained directly from the data by 
curve fitting and numerical differentiation. Then the 
fitted function and its derivative are used to 
determine the corresponding R-curve as follows: 
First, at a given point on the residual-strength curve, 
the function /, h , or j and its derivative are estimated 
by using suitable numerical methods (discussed 
later). Second, the function and its derivative are 
substituted into the appropriate version of equation 
(7), which must be solved for A c . A numerical 
solution is usually required, and in this report a 
standard program for a desktop calculator was used, 
with A c always taken as the least positive root. 
Finally A c and the fitted value of o 2 c (which equals/, 
h, or j) are substituted into equation (2) to calculate 
E'G C . The resulting pair ( E'G Ci A c ) define a point 
(E'Gr.A) on the R-curve. The process is repeated at 
additional points on the residual-strength curve. If 
desired, the R-curve can be expressed in terms of 
Kr = (E'Gr) x/1 . Also, a suitable function can be fit 
to describe Gr or Kr as an explicit function of A. 

At best this estimation method has the same 
limitation as a standard R-curve test run under load 
control; namely, for a given specimen the R-curve 
can only be measured (or estimated) to the point at 
which the crack becomes unstable. Furthermore the 
extent to which the R-curve can be estimated may 
depend on the shape of the actual R-curve, the nature 
of the residual-strength data, or both. For example, 
suppose that the actual R-curve has a pronounced 
knee (typical of high-strength, low -toughness 
materials), as shown in figure 3. There, as the initial 
crack length increases from 0.2 IT to 0.8 W , the 
location of the instability point on the R-curve moves 
only a small distance. Thus, even if one had many 
residual-strength data for such a material, having 


only normal scatter and covering the wide range 0.2 
<Xq< 0.8, the estimated R-curve plot would 
probably consist of many points grouped in a tight 
cluster. 

The best possible R-curve estimate should result 
from residual-strength data that meet the following 
criteria: 

(1) One specimen geometry factor (initial crack 
length, specimen width, or their ratio) must be 
constant, since equations (7) were derived on that 
basis. 

(2) The residual-strength data should cover a wide 
range of initial crack lengths. For constant-width, 
center-crack specimens, however, the range 
0<Xo<0.47 is adequate. 

(3) The distribution of initial-crack-length values 
should favor the shorter crack lengths. For example, 
the distribution (Xq= 0.05, 0.1, 0.2, 0.4) is better 
than (Xq =0.1 , 0.2, 0.3, 0.4). The shorter initial 
cracks help define the lower portion of the estimated 
R-curve and (because of the usual shapes of residual- 
strength curves) enhance the process of numerical 
differentiation. 

(4) The distribution of initial-crack-length values 
and the scatter in the residual-strength data should be 
such that numerical differentiation of the residual- 
strength curve can be done with confidence. 
Numerical differentiation of residual-strength data is 
discussed further in appendix A. 


Verification 

The analytical methods described herein are 
applicable to any specimen type for which a stress 
intensity analysis is available. However, residual- 
strength tests are almost always done on center-crack 
specimens. For this reason the examples to follow are 
limited to center-crack data. Experimental data sets 
containing both residual-strength and R-curve data 
from center-crack specimens are relatively few in 
number. Nevertheless enough were found to 
illustrate the application of the analysis to Case I, 
Case II, and Case III data. All calculations were 
performed by using the U.S. customary units in 
which the data were originally presented and then 
converted to SI units. 

Data for 2014-T6 Aluminum Alloy 

In reference 7 this author presented test data for 
2014-T6 aluminum alloy specimens 1.5 mm (0.06 in.) 
thick, tested at 77 K (-320° F). The specimens were 
7.5, 15, and 30 cm (3, 6, and 12 in.) wide, and initial 
crack lengths ranged from 3 mm (1/8 in.) to one- 
third of the specimen width. The average values of 
initial crack length and residual strength for the 
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longitudinal grain direction, as given in table 11(a) of 
reference 7, are used here. 

Numerical differentiation of these data was done 
as follows: Each specimen width was considered as a 
data subset, and the point (ao=0, a c = a u ) was 
treated as a common data point. A parabola was fit 
to each group of three successive data points, and the 
parabola was differentiated at the middle point. If 
there are N data points in a subset, then N-2 
deviatives can be evaluated. Also, a two-parabola 
spline function was fit to each group of four 
successive data points, as described in appendix B. 
The knot was placed midway between the second and 
third data points, and the spline function and its de- 
rivative were evaluated at the knot. In this way an 
additional N— 3 derivatives can be evaluated. Then 
equations (7a) and (2) were used to calculate points 
on the estimated R-curve, as described earlier. These 
are shown as open symbols in figure 6. 

Reference 7 also reports values of critical crack 
length. From these values the crack-extension 
resistance and the effective crack length at instability 
(failure) were computed as 


E'G r = (JcTra e ff secant (7r a e ff / W) 


(8a) 
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Figure 6. - R-curve for 2014-T6 aluminum alloy. (Data from 
ref. 7. ) 


ff eJ y=flf c +(jE"G/j/27ra^)=ao +A (8b) 

respectively. These are shown as solid symbols in 
figure 6. 

It is apparent in figure 6 that the flagged point lies 
outside the band of the remainder. This point was 
obtained after fitting a two-parabola spline function 
to four successive residual-strength data points. 
Examination of the spline coefficients shows that one 
of the fitted parabolas is concave-down. As discussed 
in appendix A this is not reasonable and hence this 
point should be ignored. The remaining points 
derived from residual-strength data by using the 
method of this report (open symbols) agree quite well 
with those calculated from critical-crack-length 
values (solid symbols). If needed, a curve could be 
fitted to these points to give an explicit equation for 
the R-curve, but this was not done here. 

Data for 2219-T87 Aluminum Alloy 

The data for 2219-T87 aluminum alloy used herein 
originally appeared in an internal report 
(Eichenberger, T. W.: Fracture Resistance Data 
Summary. Report D2-20947, Boeing Airplane Co., 
June 1962) but are also tabulated in table 1 of 
reference 8. Center -crack specimens 2.5 mm (0.10 
in.) thick were tested. The data for specimens 60 and 
120 cm (24 and 48 in.) wide are used here because 
they cover a wide range of crack lengths. 

The numerical differentiation of these data was 
done in exactly the same manner as that of the 
2014-T6 data, except that the points given by 
(ao = W/ 2, <j c = 0) were also treated as data points. 
Once again, reported values of critical crack length 
were also used in equations (8) to compute crack- 
extension resistance and effective crack extension at 
instability. Figure 7 shows that, as before, points 
derived from residual -strength data (open symbols) 
agree well with those calculated from critical-crack- 
length values. These data were also used in refer- 
ence 4, and there an explicit R-curve equation was 
needed. The simple exponential form 

E'Gr =8.07 X 10 15 A°*554 


where E'Gr is in N 2 /m 3 and A is in cm, or 
E'G r = 11.2X10 9 A 0 - 554 


where E f G R is in lb 2 /in 3 and A is in inches, provides 
a good fit to the open symbols. Residual strengths 
calculated from this equation by using conventional 
instability analysis were found in reference 4 to agree 
quite well with the original data over the entire range 
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Figure 7. - R-curve for 2219-T87 aluminum alloy. (Data from 
ref. 8. ) 


of crack lengths, with the average error being less 
than 3 percent. 

Data for AM355 Alloy Steel Sheet 

Reference 9 includes some unusual test data for 
AM355 alloy steel specimens 0.5 mm (0.02 in.) thick. 
One series of tests used specimens of various sizes 
having a constant ratio of initial crack length to 
specimen width (Case II). Another series was run 
with constant initial crack length and varying widths 
(Case III). In addition, a considerable amount of 
subcritical crack growth data is reported. 

Specimens with and without stiffeners (to prevent 
local crack buckling) were tested, but only stiffened- 
specimen data are considered here. Numerical 
differentiation of the constant-ratio (Case II) data 
(table XI, ref. 9) was done in the same manner as 
before except that the point (<x c =<r u > <*o =0) was not 
used. Then equations (7b) and (2) were used to 
calculate points on the estimated R-curve. The 
constant-crack-length (Case III) data required a 
different treatment. The residual-strength data (table 
VI, ref. 9) are plotted in figure 8. One point appears 
to be erroneously low and was ignored. The re- 
maining points do not lie along a simple, smooth 
curve, and direct numerical differentiation would not 



0 4 8 12 16 20 24 

Specimen width, W, in. 


Figure 8. - Residual strength and fitted curve for AM355 steel. 
(Data from ref. 9. ) Initial crack length, 2an, 10 cm(4 in. ) ; 
ultimate tensile strength, o u , 1660 MPa (240. 8 ksi). 
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Figure 9. - R-curve for AM355 steel. (Data from ref. 9. ) 
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be prudent. Some type of data smoothing is needed. 
The equation plotted in figure 8 is relatively simple, 
has the proper curvature and limits, and appears to 
fit the data fairly well. This equation was used to 
compute smoothed values of stress (which replaced 
the experimental values) and associated derivatives. 
These were then used in equations (7c) and (2) to 
calculate points on the estimated R-curve. Subcritical 
crack growth data are presented in figures 20 and 22 
of reference 9. These consist of four to eight points 
per specimen and were taken during the constant- 
ratio (Case II) tests. For each point the crack- 
extension resistance and the effective crack length 
were calculated by using equations (8). 

Points on the estimated R-curve, computed by 
three methods, are shown in figure 9 and are in 
generally good agreement. Those derived from the 
constant-crack-length (Case III) residual-strength 
tests lie along the high side of the scatter band, and 
this may be due to the choice of the function used to 
smooth the data. 


Concluding Remarks 

This report has presented a method for estimating 
the R-curve from residual-strength data. The method 


was verified by using various types of residual- 
strength data from the literature. Numerical 
differentiation of the residual-strength data is 
required, and the data used herein presented no 
significant problems in that respect. 

For other data not examined here the success or 
failure of the method hinges on the numerical 
differentiation process, which for residual-strength 
data is not a well-established procedure. Some 
guidelines are given in appendix A, but engineering 
judgment may still be required. If the data can be 
differentiated with confidence, the method should be 
successful. If not, the results can be misleading or 
erroneous. In that case it may be better to fit the data 
to a semiempirical fracture analysis and derive an 
equivalent R-curve by using the method of refer- 
ence 4. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, May 20, 1980, 

505-33. 
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Appendix A 

Numerical Differentiation of Residual-Strength Data 


Numerical differentiation is a mathematical 
procedure that can be done with great elaboration 
and high precision but which nevertheless can often 
be highly subjective. A set of experimental data can 
be thought of as a number of estimated values of a 
function whose analytical form is, for the most part, 
unknown. The problem is to find a simple, 
differentiable function that closely approximates the 
unknown function over at least part of the test range. 
There are two steps to the solution: The first is the 
choice of a suitable form (i.e., polynomial, 
exponential, trigonometric, etc.) for the 
approximating function; the second is the 
determination of the best numerical values for the 
coefficients of the function. The second step requires 
only established mathematical methods, but the first 
requires engineering judgment. Numerical 
differentiation is more sensitive to the choice of func- 
tional form than are interpolation or numerical 
integration. For example, consider the following 
sketch, where f(x) is the true (unknown) function and 
P(x) is the approximating function. Here the 
functional form of P(x) seems well chosen, and the 
errors in slope are small. However, a straight line will 
also fit within the bounds of /(x*)dt€ and would give 
much larger slope errors. Thus the choice of 
functional form is often the more important step in 
numerical differentiation. 

In some experiments the proper form for the 
approximating function is known in advance, either 
from first principles or from historical evidence. For 
residual -strength data there are no such guidelines. It 
seems prudent therefore to use a simple form for the 
approximating function and to check the fitted 
function for reasonableness. In this report second- 
degree polynomials (parabolas) were used almost 


exclusively. In some cases two parabolas were linked 
to form a spline function (appendix B). These simple 
functions appear to be adequate, at least for the data 
that were examined here. 

Once the chosen function has been fit to data and 
the values of its coefficients determined, the fitted 
function should be checked for reasonableness. This 
may be done by considering some basic 
characteristics of residual-strength curves. If 
specimen width is constant (Case I), the slope of the 
residual-strength curve should always be negative. 
The second derivative should be positive (i.e., the 
curve is concave-up) except as na§ approaches W and 
possibly as ao approaches 0. Also, by combining 
equations (5) and (7a) and observing figure 4 we can 
show that 


_ 1 da n 
~ a da^ > 2W ymin 

for constant -width tests. For tests where the ratio 
ao/ IF is constant (Case II) the slope should always be 
negative and the second derivative should always be 
positive. For constant-crack-length tests (Case III) 
the slope of the residual strength curve do/dW is 
always positive and the second derivative is always 
negative (concave-down). Furthermore the curve is 
always asymptotic to a finite strength value as W 
becomes large with respect to ao • Any fitted function 
that conflicts with these characteristics should be 
considered suspect, and engineering judgment will be 
required to obtain a satisfactory numerical 
differential. Graphical differentiation is a last resort, 
but because of its subjectivity it was not considered in 
this report. 
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Appendix B 

Two-Parabola Spline Function 


A spline function is a sequence of piecewise 
polynomials of degree N> 1 whose coefficients are 
such that the function itself and its first AT- 1 deriva- 
tives are continuous over its entire range. The 
junctions of the polynomials are called knots. In this 
report some of the numerical differentiation was 
accomplished by fitting two parabolas through four 
consecutive data points to form a spline function. 
The function and its first derivative were then 
evaluated at the knot, which was midway between the 
second and third points. This is done as follows. 

Let four consecutive data points be denoted as (. Xj , 
y{) where x\ <*2 <*3 <*4. We wish to fit two 
parabolas through the four points such that 


y\ = <*o+<*l*i + <*2*1 > 

yi = <*0 + <* 1*2 +<* 2*2 
,y3=0O+01*3+02*f 
.>'4=00+01*4+02*4 


(Bl) 


If the knot is located at x=Xk such that X 2 <*Ar <*3 » 
the continuity conditions at the knot are 


>/t = «O + ai*/t + a2*^ = 0o + 0l*/t + 02*jt (B2) 

y'k— a \ +2a2*Ar = 01 +202*A: (B3) 


where the prime denotes differentiation. After 
subtracting the appropriate part of equation (B2) 


from each of equations (Bl), expanding and 
rearranging terms, and substituting equation (B3), 
equations (Bl) can be written as 


y\ =^+(*1 -xk)ylc+(x 1 -**) 2 «2 + (O)02 
yi =yk + (*2 -xk)y'k + (*2 -**) 2 «2 + (O) 0 2 
y3 =yk+{x 3 -*jtM +(0)a2 +(*3 *it) 2 02 
.V4 =yk+(x 4 ~x k )y' k + (0)c* 2 +(*4 -*fc) 2 02 


which are easily solved. When differentiating 
residual-strength data, this approach has two 
advantages. First, the interpolation value y k and the 
derivative y' k are obtained directly. Second, the signs 
of the coefficients ct2 a nd 02 indicate the curvature 
of the parabolas. As discussed in appendix A, both 
oi2 and 02 should be positive for Case I or Case II 
data and negative for Case III data. If c*2 and 02 
have opposite signs, there will be an inflection at the 
knot and the derivative should be regarded with 
suspicion. If the remaining coefficients are desired, 
they are 


«i =y' k -2x k ot2 


«o =yk -x k (yk - x k a 2 ) 


0 i =y k -2x k 0 2 
0o=yk-x k (y' k -x k 0 2 ) 


li 
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